Density perturbations from both the inflaton and the curvaton 
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We consider a supersymmetric grand unified model which leads to hybrid inflation and solves the strong CP and 
fj, problems via a Peccei-Quinn symmetry, with the Peccei-Quinn field acting as a curvaton generating together with 
the inflaton the curvature perturbation. The model yields an isocurvature perturbation too of mixed correlation 
with the adiabatic one. Two choices of parameters are confronted with the Wilkinson microwave anisotropy 
probe and other cosmic microwave background radiation data. For the choice giving the best fitting, the curvaton 
contribution to the amplitude of the adiabatic perturbation must be smaller than 67% at 95% confidence level and 
the best-fit power spectra are dominated by the adiabatic inflaton contribution. This case is disfavored relative to 
the pure inflaton scale-invariant case with odds of 50 to 1. For the second choice, the adiabatic mode is dominated 
by the curvaton, but this case is strongly disfavored relative to the pure inflaton case (with odds of 10 r to 1). 
Thus, in this model, the perturbations must be dominated by the adiabatic component from the inflaton. 



1. INTRODUCTION 

The usual assumption |1I2| in inflationary cos- 
mology is that the density perturbation is gen- 
erated solely by the slowly rolling inflaton field 
and is, thus, expected to be purely adiabatic and 
Gaussian. However, lately, the alternative pos- 
sibility lj that the adiabatic density pertur- 
bations originate from the inflationary perturba- 
tions of some other "curvaton" field which is also 
light during inflation has attracted much atten- 
tion. In this case, appreciable isocurvature per- 
turbations 1 1I5| in the densities of the various 
components of the cosmic fluid as well as signifi- 
cant non-Gaussianity can arise. The main reason 
for advocating this curvaton hypothesis is that 
it makes |H] the task of constructing viable infla- 
tionary models much easier, since it liberates us 
from the restrictive requirement that the inflaton 
is responsible for the curvature perturbation. 

The standard curvaton hypothesis |3I4| . which 
insists that the total curvature perturbation orig- 
inates solely from the curvaton, can also be quite 
restrictive and not so natural. Indeed, in this 
case, one needs to impose 6 an upper bound 
on the inflationary scale in order to ensure that 
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the perturbation from the inflaton is negligible. 
This bound can be quite strong if the slow-roll pa- 
rameter e is small. Moreover, in generic models, 
one expects that all the scalars which are essen- 
tially massless during inflation contribute to the 
total curvature perturbation. So, in the presence 
of a curvaton, it is natural to assume that the 
adiabatic density perturbation is partly due to 
this field and partly to the inflaton. Finally, the 
recent measurements on the cosmic microwave 
background radiation (CMBR) by the Wilkin- 
son microwave anisotropy probe (WMAP) satel- 
lite [7| have considerably strengthened |8I9I10| the 
bound on the isocurvature perturbation. Thus, 
the viability of many curvaton models is in doubt. 

We, thus, take [TT] a more liberal and natural 
attitude allowing a significant part of the total 
curvature perturbation to originate from the in- 
flaton. We can then hopefully relax the tension 
between curvaton models and the WMAP data 
without losing the main advantage of the curva- 
ton hypothesis, which is that it facilitate the con- 
struction of viable inflationary models. 

2. THE PARTICLE PHYSICS MODEL 

In order to explore this possibility, we take [TT] 
the concrete supersymmetric (SUSY) grand uni- 
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fied theory (GUT) model of Ref. [2], which is 
based on the left-right symmetric gauge group 
G LR = SU(3) C x SU(2) L x SU(2) fl x U(1) B _ L . 
This model leads naturally to standard SUSY 
hybrid inflation |13I14| , which is the most promis- 
ing inflationary scenario. It also solves simultane- 
ously the strong CP and fj, problems via a Peccei- 
Quinn (PQ) symmetry. The PQ field, which 
breaks spontaneously the PQ symmetry, can act 
as curvaton (see Ref. |15|~) . 

The SU(2) l doublet left-handed quark and lep- 
ton superfields are denoted by % and U respec- 
tively, whereas the SU(2)# doublet antiquark and 
antilepton superfields by q\ and If respectively 
(i=l,2,3 is the family index). The electroweak 
Higgs superfield h belongs to a bidoublet (1, 2, 2)q 
representation of Glr. The breaking of G^p to 
the standard model gauge group Gsm > a t a super- 
heavy scale M ~ I0 16 GeV, is achieved through 
the superpotential 

swi = KS{r H r H - m 2 ), (i) 

where Ijj, l c H is a conjugate pair of SU(2)^ dou- 
blet left-handed Higgs superfields with B — L 
charges equal to 1,-1 respectively, and S is a 
gauge singlet left-handed superfield. The dimen- 
sionless coupling constant k and the mass param- 
eter M are made real and positive by rephasing 
the fields. The SUSY minima of the scalar poten- 
tial lie on the D-flat direction l c H — Ijj at (S) = 0, 
I (1%:) I = I Ch) I = The model also contains two 
extra gauge singlet left-handed superfields N and 
TV for solving |16| the /i problem via a PQ sym- 
metry, which also solves the strong CP problem. 
They have the following superpotential couplings: 
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Imp Imp 
where A and (3 are dimcnsionlcss coupling con- 
stants (made real and positive) and mp ~ 2.44 x 
I0 18 GeV is the reduced Planck mass. 

The model possesses three global U(l) sym- 
metries, namely an anomalous PQ symmetry 
U(I)pq, a non-anomalous R symmetry U(1).r, 
and the baryon number symmetry U(l)#. The 
PQ and R charges of the various superfields are 

PQ : q c , l c ,S,l c H , l c H (0),h,N(l), q, l,N(-l); 
R: h,l c H J c Hl N(0) iq ,q c ,l,l c ,N(l/2),S(l). (3) 



The superpotential in Eq. (JJJ leads |13ll4j nat- 
urally to the standard SUSY realization of hybrid 
inflation |17| . The inflationary path is a built-in 
classically flat valley of minima at l c H = l c H = 
and for |5| greater than a critical (instability) 
value S c = M. The constant tree-level poten- 
tial energy density k 2 M 4 on this path can cause 
inflation as well as SUSY breaking leading [H] 
to one-loop radiative corrections which provide 
a logarithmic slope along this path necessary for 
driving the system towards the vacua. 

The scalar potential Vpq for the PQ symme- 
try breaking is derived from the first term in the 
right hand side (RHS) of Eq. J2J and, after soft 
SUSY breaking mediated by minimal supergrav- 
ity, is [TT3] 
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where m 3 / 2 ~ 1 TeV is the gravitino mass and A 
is the dimensionless coefficient of the soft SUSY 
breaking term corresponding to the first term in 
the RHS of Eq. 0. Here, the phases of A, N and 
N are adjusted and \N\, \N\ are taken equal so 
that Vpq is minimized. Moreover, rotating TV on 
the real axis by a R transformation, we defined 
the canonically normalized real scalar PQ field 
(j> = 2N. For \A\ > 4, Upq has a local minimum 
at cj> — and absolute minima at 
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with f a being the axion decay constant. Substi- 
tuting this vacuum expectation value (VEV) into 
the second term in the RHS of Eq. @ , we obtain 
a fi term with 
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as desired [T^]. Note that the potential Upq in 
Eq. igj should be shifted JS] by adding to it the 
constant 
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so that it vanishes at its absolute minima. 
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3. THE PQ FIELD EVOLUTION 

In the early universe, the PQ potential can ac- 
quire sizable corrections from the SUSY breaking 
caused by the presence of a finite energy density 
[13119120] . These corrections are particularly im- 
portant during inflation and the subsequent infla- 
ton oscillations. We will ignore the A term type 
corrections [20) . To leading order, we then just 
obtain a correction Sm^ — ^H 2 (assumed to be 
positive) to the mass 2 of the curvaton, where H is 
the Hubble parameter and the dimcnsionless con- 
stant 7 can have different values during inflation 
and inflaton oscillations. We assume that Srn 2 ^ 
is (somewhat) suppressed, which can be [21] the 
case for specific (no-scale like) Kahler potentials. 
For simplicity, during inflation where the cancel- 
lation of <5m 2 can, in principle, be "naturally" ar- 
ranged to be exact (see fourth paper in Ref. [21)1. 
we take 7 = 0. During inflaton oscillations, on 
the other hand, we choose 7 = 0.1. 

After reheating, the universe is radiation dom- 
inated and, thus, H ~ l/2t < l/2t Teh = r infl /2, 
where t is the cosmic time and t le h = rri the 
time at reheating with r; n fl being the inflaton de- 
cay width. The reheat temperature T re h is [2] 

( 45 \ 3 1 
V 27r 9* J 

where g* is the effective number of massless de- 
grees of freedom (5* = 228.75 for the MSSM spec- 
trum). It should satisfy [22] the gravitino con- 
straint: T rch < 10 9 GeV. The PQ potential can 
acquire temperature corrections too both during 
the era of inflaton oscillations (from the so-called 
"new" radiation [23) emerging from the decay- 
ing inflaton) and certainly after reheating. It has 
been shown [15|. however, that these corrections 
are always subdominant and can be ignored. 

The full effective scalar potential for the PQ 
field in the early universe is given by 

Vpq = V PQ + \l 2 H 2 4> 2 + V , (9) 

whereas the full effective scalar potential V rele- 
vant for our analysis is obtained by adding to Upq 
the potential for standard SUSY hybrid inflation 
(see e.g. Ref. [2]). The evolution of (p is generally 



governed by the classical equation of motion 

4> + 3ff</> + V = 0, (10) 

where overdots and primes denote derivation with 
respect to t and the PQ field <f> respectively. This 
equation holds [24) during inflation for the mean 
field in a region of fixed size bigger than the de 
Sitter horizon provided that the quantum fluctu- 
ations of 4> do not overshadow its classical kinetic 
energy. This requires |5S] that the value 0/ of <j> 
at the end of inflation exceeds a certain value 4>q 
given by 

(11) 

where -ff; n fl is the almost constant Hubble param- 
eter during inflation. We exclude the "quantum 
regime" (i.e. 4>f < <Pq), where the calculation of 
the curvaton spectral index is not so clear. 

For <f> to receive a super-horizon spectrum of 
perturbations from inflation and thus be able to 
act as curvaton, we must make sure that this 
field is effectively massless, i.e. V" < H 2 , during 
(at least) the last 50 — 60 inflationary e-foldings. 
This requirement, which also guarantees that <j> is 
slowly rolling during the relevant part of inflation 
and emerges with vanishing velocity at the end of 
inflation, yields an upper bound on (j)f. 

The subsequent evolution of (f> during inflaton 
oscillations is given by Eq. 1)10( 1 with H = 2/3t. 
One finds |15) that, depending on the value of 4>f, 
the PQ system eventually enters into damped os- 
cillations about either the trivial (local) minimum 
of Vpq at 4> = or one of its PQ (absolute) min- 
ima at <fr — ±/ . Of course, values of 4>f leading 
to the trivial minimum must be excluded. 

The damped oscillations of 4> continue even af- 
ter reheating, where H becomes equal to l/2t, 
until 4> decays via the second coupling in the RHS 
of Eq. (J2J) into a pair of Higgsinos provided that 
their mass /x does not exceed half of the mass 
of (f> (see Ref. QS]), which is equal to 

^ (\A\ 2 - 12)* (\A\ + (W 2 -12)*)* (12) 

and is independent of A. The decay time of the 
PQ field is t$ = T7 1 , where is its decay width, 



4 



which has been found |15| to be given by 

r * - 8^4 m *- 



(13) 



Note that the coherently oscillating PQ field 
could evaporate as a result of scattering with 
particles in the thermal bath before it decays into 
Higgsinos. However, one can show ^5] that, in 
our model, this does not happen. 

4. THE CURVATURE PERTURBATION 

The perturbation 8<p acquired by (f> during in- 
flation evolves at subsequent times and, when 4> 
settles into damped quadratic oscillations about 
the PQ vacua, yields a stable perturbation in the 
energy density of this field. After the PQ field de- 
cays, this perturbation is transferred to radiation, 
which also carries a curvature perturbation from 
the inflaton. So, the total curvature perturbation 
is of mixed origin. 

The total curvature perturbation on unper- 
turbed hypersurfaces is given by 
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where p and p are the total energy density and 
pressure, and Sp the total density perturbation. 
After the curvaton decay, it becomes 5j 



C = (i - f)d + /Co, 



(15) 



where Q = 8p r /Ap r and Cc = Sp^/Sp^ are the 
partial curvature perturbations on spatial hyper- 
surfaces of constant curvature from the inflaton 
and the curvaton respectively at the curvaton de- 
cay with p r and p^ being the radiation and <fi en- 
ergy densities respectively, and 8p r and 6p,p the 
corresponding perturbations. Also, 
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evaluated at the time of the curvaton decay, and 
£ c = 2(50o /30o j where (/>o is the amplitude of the 
quadratic oscillations of (f> and 8<po the perturba- 
tion in originating from the perturbation 5(j)f 
in cf>f at the end of inflation. 



The comoving curvature perturbation 7?. ra d, for 
super-horizon scales, is given by 

(1 " f)A (lk) ^ &i + fAc ( ik) ^ ^ (17) 

where k is the comoving (present physical) wave 
number, Ha = lOOh km sec -1 Mpc -1 is the 
present Hubble parameter and &i, a c are inde- 
pendent normalized Gaussian random variables. 
Also, Ai and A c are, respectively, the amplitudes 
of — Q and — £ c at k — H , and and v c are 
the spectral tilts of the inflaton and curvaton re- 
spectively. The corresponding spectral indices are 
m = 1 + 2i^ and n c — 1 + 2v c . We do not con- 
sider running of the spectral indices, since this is 
negligible in our model. 

The amplitude Ai of the partial curvature per- 
turbation from the inflaton is ^2] 
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A(z) = (z + l)ln(l + z- 1 ) + (z-l)ln(l-z- 1 ),(19) 
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Here, Nq is the number of e-foldings suffered by 
our present horizon during inflation, xq — Sq /M 
with Sq being the value of \S\ when our present 
horizon crosses outside the inflationary horizon, 
and Xf — Sf/M with Sf being the value of\S\ at 
the end of inflation. 

The slow-roll parameters for the inflaton are [2] 

2 

zA(z) 2 , (21) 
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where z — x 2 with x — \S\/M. In the presence 
of the curvaton, the slow-roll conditions take the 
form e, |?7i| < 1, where 
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and Xf corresponds to rfc = — 1. For k -C 1, 
x/ is "infinitesimally" close to unity and we can 
put Xf = 1 in Eq. fflty. For larger values of «, 
inflation can terminate well before reaching the 
instability point at x = 1. 

Finally, k, Nq are given by 2 
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and rii = 1 — 6e + 2rji ~ 1 + 2f^ , where e and T)i are 
evaluated at the time t* when our present horizon 
scale crosses outside the inflationary horizon. 

To calculate the amplitude A c of the partial 
curvature perturbation from the curvaton, we 
take the perturbation 6(j>* = (-ff»/27r)a c acquired 
by 4> from inflation at i* (-H* is the value of H at 
i*) and follow its evolution until the end of infla- 
tion, where we find 5<j)f . To this end, we consider 
the equation of motion for <fi during inflation (see 
Eq. pOfll in the slow-roll approximation, which 
holds for the curvaton too: 

3H0 + V' = 0. (27) 

Taking a small perturbation 5<p of <f>, this gives 

3H6<p + 3H'5<j>j> + V"S4> = 0, (28) 

which using Eq. I|27|) and the Friedmann equa- 
tion, becomes 



where 
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Integration of Eq. I|29[> from t* until the end of 
inflation (at time tf) yields 

H* [ Nq 
Scpf=—a c exp (e c -T) c )dN, (31) 
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where we used the relation dN = —Hdt for the 
number of e-foldings N(k) = Nq + \n(H /k) suf- 
fered by the scale during inflation. 

For each (f>f, we construct the perturbed field 
<pf +S(j)f. We then follow the evolution of cf>f and 
4>f + 5(f> t until the time of the curvaton decay 
and evaluate Sp^/p^ at this time. The amplitude 
A c is given by 



A c Q-c — ~ 



(32) 



We have found numerically that the perturbation 
5<po in the amplitude of the oscillating curvaton at 
trj, is proportional to 5<pf. So £ c has the same spec- 
tral tilt as 5<pf, which can be found from Eq. (|31(l : 

d\nA c . . 

Vc ^ dink = ~ e ~ €c + 7/c ' ' 33 J 

where we used the relation din k — Hdt, and e, e c 
and rj c are evaluated at . The curvaton spectral 
index is n c = 1 — 2e — 2e c + 2?7 C ~ 1 + 2r; c . 

5. ISOCURVATURE PERTURBATIONS 

At reheating, gravitinos are thermally pro- 
duced and decay well after the big bang nucle- 
osynthesis (BBN) into a particle and a sparticle, 
which subsequently turns into a (stable) lightest 
sparticle (LSP) contributing to the relic abun- 
dance of cold dark matter (CDM) in the universe. 
For simplicity, we assume that there are no ther- 
mally produced LSPs. Baryons are produced via 
a primordial leptogenesis |27] occurring at re- 
heating. So, both the LSPs and the baryons orig- 
inate from reheating and, thus, inherit the partial 
curvature perturbation of the inflaton, i.e. 
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The isocurvature perturbation of the LSPs and 
the baryons is then 
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(35) 



where Clsp+b = Clsp = Cb is the partial curva- 
ture perturbation of the LSPs and the baryons 
and we used Eq. I|15fl . Here, we assume that the 
curvature perturbation in radiation (£ 7 ) practi- 
cally coincides with the total curvature perturba- 
tion. This corresponds to a negligible neutrino 
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isocurvature perturbation, which is [5] the case 
provided that leptogenesis takes place well be- 
fore the curvaton decays or dominates the energy 
density. Applying the definitions which follow 
Eq. l(T7jl . Slsp+b becomes equal to 



■if A, ( l a i + 3fAj^- 



(36) 



The model contains axions which can also con- 
tribute to CDM. They are produced at the QCD 
phase transition well after the curvaton decay and 
carry an isocurvature perturbation 



(37) 



S a = A a [ — 



where A a is its amplitude at the present horizon 
scale, v a its spectral tilt (giving the spectral in- 
dex n a = 1 + 2v a ) and a a a normalized Gaussian 
random variable independent from a,i and a c . 
The amplitude A a is given by |15j 
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where 8 is the initial misalignment angle, i.e. the 
phase of the complex PQ field during inflation, 
and (/)* is the value of </> at t*. In our case, 9 lies 
[To] in the interval [— 7r/6,7r/6] and is determined 
from the total CDM abundance Ocdm/i 2 which is 
the sum of the relic abundance [2*§| 
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of the LSPs and the relic axion abundance [50] 

1.175 

(40) 
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Here flj = Pj/p c with pj being the present en- 
ergy density of the jth species and p c the present 
critical energy density of the universe, and TOlsp 
is the LSP mass. The spectral tilt v a is evaluated 
by observing that A a depends on the scale 
only through H* /<(>*. We find 
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top 
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where the e and e c are evaluated at t* . The axion 
spectral index is n a ~ 1 + 2mp-\/2ic/<^*. 



Combining Eqs. (|36[1 and (|37|l . we find that the 
total isocurvature perturbation is |15| 
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where we used the definitions S^lsp+b = ^lsp + 
ft B and Q m = p m /p c = ^lsp+b + with p m 
being the total matter density at present. 

6. THE CMBR POWER SPECTRUM 

We will now calculate the total CMBR angular 
power spectrum Ci at the customarily used pTT] 
pivot scale fcp = 0.05 Mpc -1 . We thus define the 
amplitudes of the partial curvature perturbations 
from the inflaton (^4p,i) and the curvaton (Ap. c ), 
and the amplitude of the isocurvature perturba- 
tion in the axions (^4p, a ) at k = fcp: 



Af > : 
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The amplitude squared of the adiabatic pertur- 
bation at fcp is then given by 
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where 7?. ra d is evaluated at fcp, and the inflaton 
(Rf) and curvaton (R 2 ) contributions to R 2 are 

R 2 = (1 - f) 2 A 2 ^ and R 2 C = f 2 A 2 . c . (45) 

The curvaton fractional contribution to R is 

^ d = ^. (46) 

The amplitude squared of the isocurvature per- 
turbation at fcp is found from Eq. (|42(1 to be 



S 2 — (<5 ra d<5 ra d) — + S 2 + S 2 , 
where 

S 2 = 9 f 2 (9l^] 2 Al 



(47) 



S 2 = 9f 2 



LSP+B 



7 



5,, 



Ai 



(48) 



with 



are, respectively, the inflaton, curvaton, and 
axion contributions to this amplitude squared. 
The cross correlation between the adiabatic and 
isocurvature perturbation at fcp is 



C — (7?. ra dS, 
where 
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are the contributions from the inflaton and the 
curvaton respectively. The axions do not con- 
tribute to the cross correlation (and the ampli- 
tude of the adiabatic perturbation). Also, the 
isocurvature perturbation from the inflaton is 
fully anti-correlated with the corresponding adia- 
batic perturbation, whereas the isocurvature per- 
turbation from the curvaton is fully correlated 
with the adiabatic perturbation from it. So the 
overall correlation is mixed. It is thus useful 
to define the dimensionless cross correlation 
cos A and the entropy-to-adiabatic ratio B at fcp: 



C S 
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The total CMBR power spectrum is given, in 
the notation of Ref. |33| . by the superposition 
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These relations hold for the temperature (TT), E- 
polarization and temperature-polarization (TE) 
cross correlation angular power spectra. 

As an indicative (approximate) criterion to get 
a first rough feeling on the possible compatibil- 
ity of our model with the data, we apply the cos- 
mic microwave background explorer (COBE) con- 
straint on the TT power spectrum Cj T : 
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[(Rl + 4Sf-4a)f( ni ,£) 



+ (R 2 c +4S 2 c -4C c )f(n c ,£) 
+AS 2 J(n a ,£)], (57) 

which is accurate to about 10 — 20% in the Sachs- 
Wolfe (SW) plateau, i.e. for t < 20 (see e.g. 
Ref. Here the function f(n,£) is equal to 



T(3 - n)r(l - 1 + f ) 

2 3 -«r 2 (2 - §)r(^ + § - f) 



(58) 



.33 x 10 3 Mpc being the 



with 7/0 = 2-Hq 1 
present conformal time. 



7. NUMERICAL CALCULATIONS 

7.1. The evolution of the curvaton 

We are now ready to proceed to the numerical 
study of the evolution of the PQ field during and 
after inflation. To this end, we fix k in the su- 
perpotential 8W\ in Eq. Then, for any given 
value of Ai, we solve Eqs. (JTHJ), and 
where Xf is the solution of rjj = — 1 with rji given 
by Eq. (|2~2"j> and T Ich , which enters Eq. J5fiJ), is 
taken equal to 10 9 GeV by saturating the grav- 
itino bound 22 • We thus determine the mass 
parameter M, which is the G^ij-breaking VEV. 
Subsequently, we find the (almost constant) infla- 
tionary Hubble parameter Hi n fl = ftM 2 /V3mp. 

For any given Ai , we take a value 4>f of <fi at the 
end of inflation (at tf = 2/3H- m fi). We then solve 
the slow-roll classical equation of motion for the 
field <fi during inflation going backwards in time 
(t < t f ) by taking m 3/2 = 300 GeV, \A\ = 5 and 
a fixed A (~ 10~ 4 ) in the PQ potential. Finally, 
we put 7 = during inflation. We find that, as we 
move backwards in time, <f> increases and becomes 
infinite at a certain moment. The number of e- 
foldings elapsed from this moment until the end of 
inflation is finite providing an upper bound N max 
on the number of e-foldings which is compatible 
with <fi = (f>f at tf. 

To understand this behavior, we approximate 
the potential V by V ~ A 2 </> 6 /32mp, which holds 
for (j) — ► oo . The slow- roll equation for cf> during 
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Figure 1. The two green/lightly shaded bands in 
the Ai — (f>f plane which lead to the PQ vacua at 
t<j, for k = 3x 10" 3 , A = 10~ 4 (model A). The 
white (not shaded) areas lead to the trivial vac- 
uum and are thus excluded. The upper red/dark 
shaded area is excluded by the requirement that, 
at t*, V" < H? nfl , while the lower one corresponds 
to the quantum regime. The blue/solid line shows 
the values of Ai, <pf which approximately repro- 
duce the correct value of the CMBR large scale 
temperature anisotropy, as measured by COBE. 



inflation can then be solved analytically yielding 



4>i 



l 



4m±H h 



-At 



(59) 



where At 



t-t f 



-4m!# infl /A 2 c6 4 , 



< 0. As At -> At min = 
-> oo, which implies that the 
maximal number of e-foldings allowed for a given 
<j> f is A max ~ -iJ infl At min = 4m|iJf nfl /A 2 0/- 

We must first require that N max > Nq. The 
time t* is then found from Nq = H- m n(tf — t*). 
Furthermore, we demand that, at f*, V" < Hf n& , 
which ensures that this inequality holds for all 
times between f* and tf. This guarantees that 
<j) is effectively massless during inflation and can 
act as curvaton. It also justifies the use of the 
slow-roll approximation and ensures that the ve- 
locity of cA at the end of inflation is negligible. 



Figure 2. The two green/lightly shaded bands in 
the Ai — <pf plane which lead to the PQ vacua at 
t<f, for k = 3 x 1(T 2 , A = 1CT 4 (model B). The 
notation is the same as in Fig. ^ 



This masslessness requirement yields an upper 
bound on 0/ for every Ai and fixed k and A. 
The excluded region in the Ai — <pf plane for 
fixed k and A is depicted as a red/dark shaded 
area in the upper part of this plane. In Figs. ^ 
and [21 we show this upper red/dark shaded area 
for k = 3 x 1CT 3 , A = 1CT 4 (model A) and 
k = 3x 10~ 2 , A = 10~ 4 (model B) respectively. 
The lower red/dark shaded area corresponds to 
the quantum regime and is also excluded. 

We start from any given c6/ at tf not excluded 
by the above considerations and assume zero time 
derivative of <j) at tf . We then follow the evolution 
of 4> for t > tf by solving the equation of motion 
in Eq. gJJ with H = 2/3* for tf < t < t rch 
and H = l/2(< - t re h/4) for < rch < t < t$. 
The time of reheating i ro h is found from Eq. (JHJ) 
with = 228.75 and the curvaton decay time 
tff, from Eq. I|13fl with /3 from Eq. ©, where we 
put fi = 300 GeV. We take 7 = 0.1 after the 
end of inflation. We find that, for fixed k and A, 
there exist two bands in the Ai — <f>f plane lead- 
ing to the PQ vacua at t$. They are depicted as 
an upper and a lower green/lightly shaded band 
(see Figs, n and I2J) ■ The white (not shaded) ar- 
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eas in the A t — (fit plane lead to the false (trivial) 
minimum at (f> = and thus must be excluded. 

7.2. The calculation of C e 

For any fixed n and A, we take a grid of values of 
Ai and <f> / which span the corresponding upper or 
lower green/lightly shaded band. For each point 
on this grid, we consider <j)j and its perturbed 
value 4>j + 5(f) f and follow the evolution of both 
these fields until curvaton decay (at i^), where we 
evaluate the amplitude of Sp^f p$. The amplitude 
A c is then given by Eq. (|32Jl . For the Monte Carlo 
(MC) analysis, we use Ai and A c (rather than Ai 
and <f>f) as base parameters and limit our grid to 
Ai < 6 x 10 -5 . The indices rij and n c for each 
point on the grid are found from Eqs. I|21(l and 
(J21J) applied at x = xq, and Eqs. (23), (30) and 
(|33f) applied at i„. The fraction / in Eq. I|16[) is 
evaluated at t^. The amplitude A a of the isocur- 
vature perturbation in the axions is calculated 
from Eq. I|38|) with the initial misalignment angle 
9 evaluated, for any given total CDM abundance 
Ocdm^ 2 = Qhsph 2 + il a h 2 , by using Eqs. (f^ 
and ((40(1 with tolsp = 200 GeV. For our choice 
of parameters, the LSP relic abundance is fixed 
(~ 0.0074). The spectral index for axions is found 
from Eq. l|4T|) . 

In summary, for any fixed n and A, we take 
a grid in the variables A t and A c covering the 
upper or lower green/lightly shaded band. For 
any fish 2 and Vl a h 2 , we then calculate the am- 
plitudes squared of the adiabatic and isocurvature 
perturbations from Eqs. (@3J, (@2J and Eqs. l(4"T|) . 
1|48JI respectively, the cross correlation amplitude 
from Eqs. igHl), (EH and the total CMBR TT 
and TE power spectra via Eqs. fB^ — lf55| . The 
curvaton fractional contribution to the adiabatic 
amplitude F^ d , the dimensionless cross correla- 
tion cos A and the entropy-to-adiabatic ratio B 
are found from Eqs. 1(46(1 and (51) • ^ n computing 
M, A c and A a , we fix the present Hubble param- 
eter to the Hubble space telescope (HST) central 
value Hq — 72 km sec -1 Mpc -1 which has 
an impact less than 3% on our results. Clearly, 
we do allow Hq to vary in the MC analysis. 

Before deploying the full MC machinery to de- 
rive constraints on the parameters, it is instruc- 
tive to obtain a first rough idea about the vi- 



ability of our model by using the approximate 
expressions in Eqs. (|57|l and 1)58(1 for the temper- 
ature SW plateau and requiring that the COBE 
constraint in Eq. (|5(j|) is fulfilled. To this end, 
we take VL B h 2 = 0.0224 and fi C DM^ 2 = 0.1126, 
which are the best-fit values from WMAP [7] . We 
find that the COBE constraint cannot be satis- 
fied in the lower green/lightly shaded band in the 
Ai — cj)f plane for any k and A. The reason is 
that the relatively low values of (/)* in this band 
combined with the sizable relic abundance of the 
axions leads to an unacceptably large axion con- 
tribution to the RHS of Eq. 1(57(1 . In the upper 
green/lightly shaded band in the Ai — cf>f plane, 
on the contrary, the COBE constraint is easily 
satisfied. The resulting solution is depicted by a 
blue/solid line (see Figs, d and EJ. 

7.3. The MC analysis 

We proceed to the full MC analysis confronting 
the predictions of our model with the CMBR 
data. We use a version of the Markov chain MC 
package COSMOMC as described in Ref. |37| . 
The adiabatic and isocurvature CMBR transfer 
functions are computed in two successive calls 
similarly to Ref. jHE]- For fixed k and A, the ini- 
tial conditions are fully specified by Ai and A c . 
The MC sampling takes as free parameters the 
Ai and A Cl the present baryon and axion abun- 
dances lob = f^s/i 2 and w a = VL a h 2 , the present 
dimensionless Hubble parameter h, and the red- 
shift z r at which the reionization fraction is a 
half. All other derived quantities are computed 
from the above parameter set. In particular, the 
total CDM abundance is wcdm = ^cdm'i 2 = 
^lsp + w a , where w L sp = ^lsp^ 2 ~ 0.0074 = 
constant. Also, since we take flat cosmologies 
only, the cosmological constant energy density 
J7a (in units of p c ) is a derived parameter, i.e. 
il\ = 1 — (wcdm +ws)//i 2 . The gravitational 
waves are negligible. In summary, for fixed n and 
A, our parameter space is six dimensional: 

{u> B ,u a ,h, z r ,Ai, A c }. (60) 

We compare the predicted CMBR TT and TE 
power spectra to the WMAP first-year data |Zj 
with the routine for computing the likelihood sup- 
plied by the WMAP team At a smaller angu- 
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Figure 3. ID marginalized posterior distribu- 
tion for model A (black solid line) and model 
B (cyan/light gray solid line) with only the up- 
per green/lightly shaded bands included. The 
red/medium gray line is for the pure inflaton case 
with ni = 1, plotted here for comparison. Model 
A is quite close to the pure inflaton case. Model B 
displays a preference for non-zero curvaton con- 
tribution. Its quality of fit is, however, poorer. 
We also display as dotted smoothed curves the 
values of the mean posterior. All the curves are 
normalized at their peak value. 



lar scale, we add the cosmic background imager 
(CBI) and the decorrelated arcminute cos- 
mology bolometer array receiver (ACBAR) 
band powers as well. We use flat top-hat priors 
on the base parameters 

uj B ,uj a ,z r ,Ai,A c . (61) 

The limits of the top-hat priors do not matter 
for parameter estimation, as long as the poste- 
rior density is negligible near the limits. How- 
ever, the prior range of the accessible parame- 
ter space plays an important role in computing 
the Bayes factor for model testing. We limit the 
maximum range of h by imposing a top-hat prior 



Figure 4. As in Fig. but for (most of) the 
derived parameters (the pure inflaton model is 
not included). The adiabatic amplitude in model 
B is dominated by the curvaton (F c ad ps 0.8). 



0.40 < h < 1.00 and we use the HST result [33] 

L^ocexp^-^) 2 , (62) 

where ho = 0.72 and a = 0.08, by multiplying 
the likelihood function for the CMBR data by 
the above Gaussian likelihood. By trying differ- 
ent priors, we find JT] that the broad lines of 
the constraints for the PQ model are robust with 
respect to the choice of non-informative priors. 

We will examine in some detail the models 
A and B, which are representative cases of the 
two possible behavior patterns of our PQ model. 
We first consider the upper green/lightly shaded 
band of these models. The 1-dimensional (ID) 
marginalized posterior distributions for the base 
and most of the derived parameters are plotted, 
respectively, in Figs. |3] and 

In the upper band of model A, the power spec- 
tra are dominated by the adiabatic inflaton con- 
tribution. The quality of the fit, as expressed by 
the maximum likelihood, is slightly better than 
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for the pure inflaton case with rij = 1, which is not 
surprising since the curvaton contribution plays 
a modest role. The CMBR data yield the bound 
A c < 43.2 x 1(T 5 at 95% confidence level. The 
total TT power on large scales is slightly larger 
than the pure adiabatic part, which increases the 
height of the SW plateau relative to the height of 
the first adiabatic peak. This mimics the impact 
of a larger z r , and explains why model A prefers 
a later reionization than the pure inflaton case. 

The upper band of model B exhibits a prefer- 
ence for a non- vanishing curvaton amplitude with 
very high significance (A c — 62.8 ± 10.4). The 
overall quality of the fit is though worse than for 
model A (upper band) because this model does 
not reproduce with enough precision the shape of 
the first acoustic peak. Furthermore, this model 
shows a preference for a rather high baryon abun- 
dance (lob ~ 0.026), which is in strong tension 
with the value indicated by BBN together with 
observations of the light elements abundance, 
which typically yields u B ~ 0.020 ± 0.002 |4*2|. 

For the lower band of model A, the quality of 
the best fit is poor, whereas the lower band of 
model B is incapable of producing a spectrum in 
reasonable agreement with the data. Thus, the 
lower bands of these models are excluded. 

So far we have derived parameter constraints 
for models A and B from the data. The next 
question is to compare the upper bands of these 
two models with the pure inflaton scale-invariant 
model and decide which of these three models is 
most favored by data. Model comparison is a dif- 
ferent issue than parameter extraction and it can 
very well be that it arrives at a different conclu- 
sion. Indeed, it can be that the estimated value 
of a parameter under a model Mi is far from 
the null value predicted by model M2, but Mi 
is disfavored against M2 by Bayesian model test- 
ing. This is the case for A c in model B (upper 
band) compared to the pure inflaton model. 

To apply Bayesian model testing, we compute 
the Bayes factor B for models A and B (upper 
bands) comparing each of them to the pure in- 
flaton model with = 1. We find ^ that 
h\B = —3.9 and —17.1 in the two cases respec- 
tively. Thus, model A (upper band) is disfavored 
against the pure inflaton model with odds of 50 



to 1, while model B (upper band) with odds of 
10 7 to 1. So, model B (upper band), which pre- 
dicts a non-zero curvaton contribution to the adi- 
abatic perturbation, is strongly disfavored. On 
the contrary, model A (upper band), which does 
not require such a contribution, is only mildly 
disfavored. However, in view of the fact that it 
addresses many issues lying outside the scope of 
the pure inflaton model, such as the strong CP 
and /i problems, the generation of the BAU and 
the nature of CDM, we can by no means reject it. 

8. CONCLUSIONS 

We studied a SUSY GUT model solving the 
strong CP and /1 problems via a PQ symmetry 
and leading to standard hybrid inflation. The 
PQ field acts as curvaton contributing to the cur- 
vature perturbation together with the inflaton. 
The model predicts isocurvature perturbation too 
of mixed correlation with the adiabatic one. We 
confronted the predictions of the model with the 
CMBR data from WMAP, CBI and ACBAR and 
found that a non-zero curvaton contribution to 
the adiabatic perturbation is not favored. 
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